Possible Exciton Bose Liquid in a hard-core boson ring model 
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We present a Quantum Monte Carlo study ol a hardcore boson model with ring-only exchanges 
on a square lattice, where a K\ term acts on 1x1 plaquettes and a K2 term acts on 1x2 and 2x1 
plaquettes. At half-fiUing, the phase diagram reveals charge density wave for small K2, valence bond 
solid for intermediate K2, and possibly for large A'2 the novel Exciton Bose Liquid (EBL) phase 
first proposed by Paramekanti et al [Phys. Rev. B 66, 054526 (2002)]. Away from half-fiUing, the 
EBL phase is present aheady for intermediate K2 and remains stable for a range of densities below 
1/2 before phase separation sets in at lower densities. 

PACS numbers: 71.10.Pm, 75.10.Jm, 75.40.Mg 



The search for exotic quantum phases has gained a 
very wide audience in recent years. There has been much 
interest in "Bose-metal"-type phases in boson and spin 
systems, and more broadly non-Fermi liquids in electron 
systems. The interest has also spread to the High Energy 
community where "holographic metals" have been much 
investigated. To date, there is still no microscopic model 
where a Bose-metal phase can be controUably demon- 
strated. The novel Exciton Bose Liquid (EBL) proposed 
a number of years ago by Paramekanti et a?[l| provides 
an interesting critical bosonic phase which shares many 
characteristics normally associated with electrons in a 
metal [2] and which can be viewed as a Bose-metal in 
a restricted sense. Specifically, it is an example of a 
2D quantum system with microscopic bosonic degrees of 
freedom that nevertheless has surfaces of low-energy ex- 
citations and metal-like properties such as large specific 
heat and thermal conductivity. The striking proposal 
of Paramekanti et aZ[l| stimulated a number of works 
seeking to establish the stability of the EBL phase in 
hard-core boson models with ring exchange interactions 
on the square lattice [3|-|6[. However, these studies found 
that the EBL is not realized in these models. Instead, 
such ring interactions favor a (tt, tt) charge density wave 
(CDW) in the half-filled case, while away from half-filling 
they induce strong tendencies to phase separation. 

Here we propose a ring-only hard-core boson model on 
the square lattice with additional ring exchanges as a can- 
didate model for realizing the EBL. We define plaquette 
exchange operators 



p: 
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where by annihilates a boson on a site r, and x, y arc the 
unit vectors on the square lattice. The Hamiltonian is 



H=-K,Y. Pr" - K2 J2 (Pr' + Pr'] 



(2) 



Throughout we set Ki = 1. The original ring model 
proposed in Ref. [ij and studied numerically in Refs. Isj-la 
corresponds to K2 = 0. The extra ring terms frustrate 



the CDW order preferred by the Ki term and may allow 
the EBL phase to be realized for large K2 at half-filling. 
The present K1-K2 model has the same lattice symme- 
tries and boson number conservation on each row and 
column as the original model. From the outset, we de- 
fine our Hilbert space as the sector with equal number 
of bosons in each row and column. Note that we are al- 
lowed to make such a sector restriction consistent with 
the conservation laws of the Hamiltonian, and it is in 
this setting that we look for a uniform "liquid" phase of 
bosons that does not break any lattice symmetries. The 
original reference [l| envisioned the EBL that would be 
stable also in the presence of boson hopping. While we do 
not think that our K1-K2 model will achieve this, realiz- 
ing the EBL even in the restricted sense is already very 
interesting as it is expected to have many unusual prop- 
erties described earlier, which can then be confronted in 
a concrete 2D system. 

For simulation on large two-dimensional lattices, the 
Stochastic Series Expansion (SSE) has been the method 
of choice. Although the model docs not have a sign 
problem, implementation issues have so far precluded the 
study of ring-only models using SSE l3|. We instead ap- 
ply the Green's Function Monte Carlo (GFMC) method 
using the bias-controlled technique of Ref. [7| for an ex- 
ploratory smaller scale study. To obtain an initial wave- 
function for the GFMC projection, we first perform a 
Variational Monte Carlo study using EBL-inspired wave- 
functions. Specifically, consider a quantum rotor model 
which is a soft-core version of the iiTi-only model: 
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with canonically conjugate phase 0r and boson number 
rir operators. A harmonic "spin wave" theory where we 
expand the cosines motivates the following wavefunction, 
which we now restrict to the Hilbert space of the hard- 



core boson model 
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In the Xi-only spin wave theory, W^(q) — \J\J jK. For 
the K\-Ki case, VF(q) becomes a q-dependent function 
with two variational parameters. The formal EBL wave- 
function is quite interesting by itself — e.g., it can pro- 
duce both liquid and solid phases [9j. Here it is only 
important that the optimal trial wavefunctions provide 
good starting points for the GFMC projection. We tested 
our GFMC setup against exact diagonalization calcula- 
tions for a 6x6 lattice and obtained complete agreement 
for all physical quantities considered in this work. 

To characterize the system, we measure the density 
structure factor 

5(fe, 9«) = ;^ E ^"'■^"""""'^ ("rrv - f?) , (6) 

r,r' 

and the plaquette structure factor 

r,r' 

where {Pl^Y equals 1 if the 1x1 plaquette is hoppable 
and otherwise. This operator is diagonal in the boson 
number basis and is easier to implement in the GFMC. 
While quantitatively different from the off-diagonal P^^ 
plaquette structure factor used in Ref. y, both give qual- 
itatively the same access to bond-solid-type phases. 

We now summarize our main results at density 1/2. 
For K2 = 0, the GFMC study confirms the staggered 
CDW phase anticipated by Paramekanti et a/[l[ and ob- 
served by Sandvik et aZ[3[ in the J-K model at Kj J = 64, 
where J is the nearest neighbor hopping strength. We 
find that as K^ increases, the ground state changes to a 
Valence Bond Sohd (VBS) at K2 ^ 0.4. Interestingly, 
Sandvik et a/[3| observed a striped bond-plaquette order 
in the J-K model for 8 < KjJ < 14. (The two findings 
may be connected if, in a perturbative picture in the in- 
sulator, J induces significant K2 ring exchanges.) As 
K2 increases further, our GFMC study suggests that the 
bond order vanishes for K2 > 4, thus allowing the possi- 
bility for a realization of the EBL phase at half-filling. 

On the other hand, for a window of density below 
1/2, we already find a clear EBL phase for a wide K2 
range including also intermediate K2 values. At still 
lower density the system has strong tendencies to phase- 
separate, which occurs since ring exchanges are operative 
only when bosons are sufficiently close p, |6| . Below, we 
present evidence for these conclusions. 







K2/K1 

FIG. 1: [Color online] Structure factors versus K2 for periodic 
lattices with length L — 6,8, 10 and 12, where open symbols 
are for the density structure factor 5'(7r,7r), and solid symbols 
are for the plaquette structure factor P{it, 0). Note that linear 
scale is used for K2 = to 1 in steps of 0.1, and log scale for 
K2 = 1 to 10 in steps of 1. 
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FIG. 2: [Color online] Top: Finite size scaling of the stag- 
gered CDW order parameter S{tv,tv)/L at selected K2- Bot- 
tom: Finite size scaling of the bond-solid order parameter 
P(7r,0)/L^ at selected A'2. 



Figure [T] shows the density structure factor S'(7r, tt) and 
the plaquette structure factor P(7r, 0) plotted against K2 
for lattice sizes ranging from i = 6 to 12. At small K2, 
S{tt, tt) increases strongly with L while the size depen- 
dence becomes weaker at intermediate K2. This coin- 
cides with a strengthening size dependence of P(7r,0). 
Figure [2] shows the finite size scaling of the respective or- 
der parameters. For K2 < 0.4, the CDW order is present 
while the plaquette order is absent. For K2 > 0.4, the 
situation is reversed. Thus we conclude that the system 
undergoes a transition from the (tt, tt) CDW at small K2 
to the striped (tt, 0) [or (0, tt)] plaquette order at inter- 
mediate K2, where the latter is consistent with either a 
columnar VBS or a plaquette state. 

To get a complete picture, we measure the structure 
factors over the entire Brillouin zone, paying attention 
also to long-wavelength behavior near the lines qx = 
and Qy = 0. The CDW is an incompressible Mott insu- 
lator and is expected to have analytic S{qx,qy)- This is 
indeed what we observe [9[ and is in contrast with the 
compressible EBL behavior described below. The Mott- 
likc dependence of S{qx,qy) at long wavelengths persists 
also after the (tt, tt) CDW Bragg peak disappears for 
A'2 > 0.4, thus indicating that the system is still not 
a liquid but a different solid (here bond-solid as deter- 
mined by the corresponding order parameters). 

For the plaquette structure factor P{qj:,qy), we find 
that (tt, 0) and (0, tt) are the only sharp peaks. We also 
measure bond energy structure factor (not shown) and 
again find only (tt, 0) or (0, tt) order. This is more con- 
sistent with the columnar VBS state, which is similar to 
the phase in the J-K model[3|, Ij]. 

At larger K2 > 4:, the plaquette structure factor is 
essentially insensitive to lattice size, and we conclude, 
within the limitation of our small lattice study, that the 
(tt, 0) plaquette Bragg peak is absent. This suggests that 
the bond-solid ordering exists only at intermediate K2 
and opens the possibility of the EBL phase for large K2- 

Our present evidence for the EBL comes from the ab- 
sence of any charge or bond order. Specifically, we show 
the structure factors S(qx,qy) and P{qx,qy) over the en- 
tire Brillouin zone. Top panel in Fig. [3] shows the density 
structure factor at K2 = 5 and clear absence of any CDW 
ordering. Middle panel in Fig. [3] shows the plaquette 
structure factor for the same system, which again does 
not show bond or plaquette ordering. Despite the poten- 
tial instability hinted by the P(7r,0) and P(0,7r) cusps, 
the nearly size independence of the plaquette structure 
factor illustrated along the cut qx = i^ shown in the bot- 
tom panel of Fig. [3] gives us some confidence that there 
is no bond-solid order, so the EBL phase may indeed be 
realized in the large K2 regime at half-filling. 

Furthermore, we highlight the presence of the long 
wavelength EBL signature in the density structure factor 
near the lines g^ = and g^ = 0. The EBL theory [Ij 
predicts a V-shaped singularity [10[ for small qx across a 
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FIG. 3: [Color online] Top: Density structure factor S{qx,qy) 
for a 12x12 lattice at K2 = 5. Middle: Plaquette structure 
factor P{qx,qy) for the same system. Bottom: Cut of the 
plaquette structure factor at qx = n for lattice sizes L = 
6, 8, 10 and 12. 



cut at a fixed qy, 

<5'EBL(g2: ^ 0,qy) = c{qy)\qx\ 



(8) 



The slope c{qy) is non-zero except in the limit of small 
qy where c{qy) ^ \qy\. The EBL behavior is qualitatively 
different from the incompressible non-singular behavior 
observed in the CDW and VBS regimes. Our present 
data in the A'2 > 4 regime shows large and healthy slopes 
c{qy) consistent with the EBL. However, we caution that 
there is some weak downwards renormalization of the 
slopes upon increasing L and so far we cannot rule out 
the possibility of the EBL behavior disappearing at much 
larger lattice sizes, which would then hint that there is 



bond-solid order with very small order parameter or some 
other instability. 

We extend the search for the EBL phase to densities 
n < 1/2. Previous studies [^-Q of the J-K model have 
found that ring interactions induce strong tendencies to 
phase separation. For large K, uniform superfluid ex- 
ists only for n > 0.4, while the bosons phase separate 
at lower densities. The more extended K2 ring inter- 
actions can somewhat offset this tendency and stabilize 
a uniform liquid over a larger range of densities below 
half-filling. Here we discuss our results for two values 
K2 = 0.5 and 1.0 on a 12x12 lattice. We consider the 
model with densities varying from 2 up to 5 bosons per 
row and column (corresponding to densities n = 1/6 up 
to 5/12). Note that despite the condition of equal num- 
ber of bosons in each row and column, the system can 
still phase-separate where bosons clump, e.g., along the 
x + y direction in our samples. (Of course, in a finite sys- 
tem, this condition frustrates somewhat the tendencies 
to phase separation, but we hope that already with our 
sizes we can observe the local energetic preferences.) For 
2 and 3 bosons per row at K2 =0.5 and 2 bosons per row 
at A'2 = 1.0, we indeed see such signatures of the phase 
separation in real space and also in momentum space, 
where the density structure factor develops strong peaks 
at the smallest wavevectors 27r/L [9|. On the other hand, 
for the higher boson numbers per row, we do not see any 
signs of the phase separation or some order. The den- 
sity structure factor near q^^ = or q.y = lines shows 
clear signatures of the EBL as discussed after Eq. (|8]) 
with robust and essentially size- independent slopes [9J, 
i.e., we do not see the weak downward renormalizations 
of the slopes that constitute one nagging worry about 
our EBL for large K2 at half-filling discussed earlier. Al- 
though our study is performed on small lattices only, it 
strongly suggests that the EBL phase is realized for den- 
sities 1/3 < p < 1/2 for intermediate K2 values, while 
the phase separation sets in at lower density. 

We can rationalize our findings using the EBL theory 
PJ, Illl4l4l | . In a longer paper [9[ , we will present an in- 
teresting perspective on the EBL using slave particle ap- 
proach with two partons propagating one-dimensionally 
(ID) along one or the other lattice directions, which in 
particular allows application of quasi-lD thinking to this 
problem. The EBL fixed point theory is characterized by 
a "stiffness" and the stability requires this stiffness to be 
sufficiently large. The condition is particularly stringent 
at half-filling because of allowed Umklapp interactions. 
When the system is unstable, it is driven to either CDW 
or plaquette solid depending on the sign of some effective 
interactions. A crude scenario of what may be happen- 
ing in our K1-K2 model is that for small to intermediate 
K2 the system is unstable to Umklapps. For small K2, 
the effective nearest-neighbor repulsion (due to avoidance 
preferred by the Ki ring exchanges) drives the system 



to the CDW state. Increasing K2 effectively introduces 
a competing second-neighbor repulsion and switches the 
instability to the plaquette solid (our measurements sug- 
gest that a columnar VBS is realized instead, but the 
two bond-solids are usually closely related). Overall, in- 
creasing K2 makes the bare stiffness large, as we observe 
from monitoring the density structure factor near the 
lines ga; = or ^2^ = 0, and for large enough K2 the 
EBL becomes stable even at half-filling. When we step 
away from half-filling, the Umklapp terms are no longer 
allowed, and we can get a stable EBL already for inter- 
mediate K2- 

To conclude, our study of the K1-K2 model shows the 
transition from the CDW to VBS order at intermediate 
K2 and tantalizingly suggests the EBL phase at large K2. 
Further studies of larger sizes and confrontation with the 
theory of this unusual phase will certainly be very inter- 
esting. Even if the K1-K2 model has only solid phases, 
we have come much closer to the EBL and small addi- 
tional terms should be able to further stabilize it. The 
EBL phase is more robust away from half-filling and is 
realized in the K1-K2 model already for intermediate K2- 
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